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Abstract 

We systematically study relativistic two-species fermions with tunable attractive in- 
teractions and number-density asymmetry at zero temperature. In general, a Bardeen- 
Cooper-Schrieffer (BCS) - Bose-Einstein Condensation (BEC) - relativistic BEC (RBEC) 
crossover is observed. A generalized BCS ground state and its stability are analyzed. 
The homogeneous superfluid phase can become unstable and we consider phase sep- 
aration in real space for neutral systems. In the nonrelativistic limit, our results are 
consistent with well-known results. In addition to a BCS-BEC crossover similar to that 
in the nonrelativistic case, in the strongly attractive regime gapless excitations of an- 
tifermions and a RBEC state are observed. We address how different phases respond to 
number density asymmetry and present predictive phase diagrams of BCS-BEC-RBEC 
crossover. 
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1. Introduction 

It is well-known that in nonrelativistic systems, by tuning attractive interactions 
between fermions, a BCS superfluid of degenerate fermions can be smoothly connected 
to a Bose-Einstein condensate (BEC) of fermion pairs. For nonrelativistic Fermi gases, 
it was shown 1 1, 2] that the BCS ground state captures important features in both limits. 
Study of intermediate phases in between these two limits is called BCS-BEC crossover, 
which has been intensely studied in high-temperature superconductors 10,3] and ultra- 
cold Fermi gases | s', '6', v^] . In ultra-cold Fermi gases, it was demonstrated that number 
density asymmetry [8, 9] or heteronuclear pairing ifioll can also be realized. 

Recently, the idea of BCS-BEC crossover has been applied to relativistic systems 
of fermions ll lLll2lll3lll4lll5L ll6l]. This work is motivated by two observations. The 
first comes from color superconductivity of quarks, which has drawn intense atten- 
tion in high-energy physics UtIITsIi . At very high baryon density, the ground state of 
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quantum chromodynamics (QCD) is believed to be color superconducting states such 
as the color-flavor locking (CFL) state or the two-flavor superconducting (2SC) state 
1 18]. These states have diquark pairs, which are similar to Cooper pairs of electrons 
in conventional superconductors. At low baryon density, the system becomes hadronic 
matters fl^ I20I] . It was shown flT] that as baryon density decreases, diquark pairs 
with small pair size comparable to quark-quark distance may emerge. The conden- 
sate of these tightly bound diquark pairs are similar to the BEC of Cooper pairs in the 
BCS-BEC crossover of nonrelativistic fermionic systems. Therefore there may exist a 
crossover from a color superconducting state with weakly bound quark pairs (of large 
size) at high baryon density to a (relativistic) BEC state of diquarks at low baryon den- 



sity before the system becomes hadronic matters 0121 12 IL 11511 . Perturbative methods, 
which can only be applied at high baryon densites, are not appropriate for studying this 
crossover 

The second observation comes from the fact that a relativistic BEC (RBEC) state 
may appear even in a nonrelativistic Fermi system if the attractive interaction is strong 
enough il4ll . In an RBEC phase, there is an appreciable density of antifermions. Typ- 
ically, when the s-wave scattering length a of a nonrelativistic system approaches the 
limit a ^ {^/2m)^^, relativistic effects will become significant since the binding en- 
ergy of a fermion pair is ^ f? j (2ma}). Here m is the fermion mass. In addition to 
these observations, a crossover in QCD similar to a BCS-BEC crossover of nonrela- 
tivistic fermions may also be possible in a pion superfluid or in the interior of compact 
stars imilailll- 

We study possible phases of a two-species relativistic Fermi gas with presumably 
tunable attraction. Our approach is inspired by BCS-Leggett theory of nonrelativistic 
Fermi gases. Number density asymmetry (or polarization), which has been intensely 
studied in ultra-cold Fermi gases, is included in our calculations. This is important 
because many interesting phases such as polarized superfluid, phase separation in real 
space, Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) state |26], may be observed in non- 
relativistic neutral Fermi gases [27.1. To compare with those results, we only study 
neutral relativistic Fermi gases. 

Past work on relativistic BCS-BEC crossover introduced models including one 
species of fermions |15] as well as models including two-species of fermions 128 , 
la, 29ll . Considerations of two (or more) species of fermions are natural since pair- 
ing among quarks with different color or flavor indices may exist at moderate baryon 
densities. Ref. [28J implemented a variational approach with a coherent state. Missing 
from this work is an analysis of the stabiUty of superfluid phases. Ref. [16] studied a 
model with auxiliary boson fields. Conversion between fermions and bosons induces 
indirect fermion-fermion attraction. Although such a boson-fermion model is possible 
for ultra-cold fermions across a Feshbach resonance, it may not be necessary for pro- 
viding a mechanism for diquark pairing. Here we consider effective attraction between 
fermions without any auxiliary boson field. 

We will present phase diagrams as a function of polarization and the strength of 
the attractive interaction. Following Ref. fs^], we analyze the mechanical stability of 
homogeneous polarized superfluids. When they are unstable, one expects a phase in 
which a superfluid and a normal Fermi gas can be found in separated spatial regimes. 
This is an example of phase separation (or mixed states). This inhomogeneous phase 
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has been discussed in neutral fermionic systems such as uhra-cold fermions fsT] and 
quark matter in neutron stars 125 1 . However, phase separation may not be stable for 
charged fermions 13211 . In charged Fermi gases with number density asymetry, there is 
another instability, magnetic instability (or chromo-magnetic instability in color super- 
conductivity), which is associated with the existence of an imaginary Meissner mass 
of gauge bosons 13311 . In the mechanically unstable regimes of homogeneous super- 
fluids, there are other possible phases in additional to phase separation. For example, 
the LOFF state [26.1 is also possible in the weakly attractive regime, but here we do 
not consider this state since we are mainly interested in the strongly attractive regime. 
A third possibility are chiral condensates with (ipip) ^ 0, which become dominant 
when baryon density is low enough. A phase of competition and/or coexistence be- 
tween (R)BEC states of i/'V'-condensates and ?W;-condensates may exist before the 



chiral condensate phases completely dominate 1211 119L 12011 . We do not address these 
more complicated situations here. Collective excitations beyond mean field theory are 
also interesting and some discussions can be found in |29]. 

This paper is organized as follows. In Sec.|2]we present the basic theoretical frame- 
work of our mean field theory for two-species Fermi gases with fermion-fermion inter- 
actions. We implement both a Nambu-Gorkov formalism and a Bogoliubov canonical 
transformation and derive equations of state. In Sec. [3] we address the important sta- 
bility criterion by considering properties of the number susceptibility matrix. Phase 
separation, which is possible in the unstable regimes of homogeneous superfluids, is 
also discussed. In Sec.|4l we present our numerical results and discuss their implica- 
tions. 



2. Theoretical Framework 

2.1. Nambu-Gorkov Formalism 

We consider the following Lagrangian 

L= ^ 4'Sl^c>f,~ma + ^i„-1°)^„+Li, (1) 

'y=^,^ 

where ipa, are Dirac four-component spinors with fermion masses to^, and Lj 
denotes the attractive interactions between the two species of fermions in the = 
O"*" channel ll34i [35ll . The fermions are assumed to be neutral. The metric is 5^'' = 
diag(l, —1,-1, —1). We use the spin indices cr =|, | to denote the two species. In this 
paper we approximate = — m. The interaction Lj is of the form 

Li = -.g(7/.f «75^i)(V5ii75CV;f ), (2) 

where g is the attractive coupling constant and C ~ ij'^j'^ is the charge conjugation 
matrix. We introduce the gap function A = —g{tp'^Ci'jc,^i), which corresponds to the 
order parameter. From mean field approximation the Lagrangian becomes 

Lmf = '/'-(^7^9^ -m + Ma7°)V'<T + A*{^JCij5^i) + A(Vi;Z75Ci^f ). (3) 
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To define finite-temperature Heisenberg operators, we take t — it and x — (t, x). 
Similar to the BCS th eory of superconductivity, here we introduce a Nambu-Gorkov 
spinor following Ref . lll4ll : 



We also introduce a Green's function matrix 

S{X,X') = -{Tr[^ix)^{x')]). 

Here denotes T-order. From equations of motion one has 

D{x)S{x,x') = 18x8(5(2; - x'), 



(4) 



(5) 



(6) 



where 6{x — x') — 5{t — t')(5(x — x'), l8x8 is the eight-dimensional identity matrix, 
and D{x) is a differential operator matrix. Its expression can be found in Appendix lAl 
In momentum space, one has 



D{K)S{K) = l8x8, 



(7) 



Here K — {iutn, k) is a four-momentum at finite temperature with uin = (2ri + 1)tt/(3 
denoting the fermion Matsubara frequency. Here we set fc^ = 1 and therefore (3 = 
1/T. We present our formalisms at general temperatures but only consider results at 
T = 0. 

The inverses of D{K) corresponds to interacting fermion propagator matrix 



where the matrix elements are 

Similar expressions for G'|(i^, — and F||(iir, — can be found. The 
fermion self-energy is defined as 



(8) 

(9) 
(10) 



(II) 



A similar expression for I]|(i4r, — /i^) can be found. The energy dispersions of relativis- 
tic fermions and antifermions are = ± fi^r, where = \/k-^ + m^. Here "— " 
is for fermions and "+" is for antifermions. The Fermi energy Ep = y^kp + mi? and 
the Fermi momentum kp of an unpolarized non-interacting Fermi gas with equal total 
number density are chosen as units of energy and momentum. Here n = fc|,/(37r^) 
and n denotes total number density. 
We introduce two energy projectors 



A±(k) 



1± 



7^(7 • k + m) 



(12) 
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By using these energy projectors, the non-interacting femiion propagator matrices can 
be written as 

^ A+(k) , A_(k) 



IU!„ 



7 



(13) 



and a similar expression for Gqi {K, One can see that propagator of each species 
has contributions from both fermions and antifermions. The interacting fermion prop- 
agator matrices can be expressed as 



,(-)2 



A+(k) 



ka 



,(+)2 



,(-)2 



A+(k) 



+ - 



A_(k) 



(+)2 
k 



^kfS 

A_(k) 



tU), 



+ Et 



ka 



IUJ„ 



^k/S 



(14) 



A similar expression can be found for Gi{K, ^i, — /ti|). In the above expressions we 



introduced energy spectra for fermions, E,.^ ^ 
^ Here we define EZ - 



Eka,l3 



Ek ± Cfefc- 



2 



± 



,(±)2 



k k 

the case of Fermi gases with only one species, here both energy spectra of fermions 
and antifermions may cross zero when there is a population asymmetry in the two- 
species. As a consequence, there may exist gapless excitations for both fermions and 
antifermions. We define fi — (fi^ + fii)/2 and h = (/if — IJ-i)/2. We choose the 
spin-up species as the majority so /i > 0. There are zeros in energy spectra of fermions 
and antifermions. Zeros of the energy spectrum of fermions, if exist, are at fc^ 2 = 
max(0, (jU ± \/h? — A^)^ — m?). Zeros of the energy spectrum of antifermions, if 
exist, are at fc| 4 = max(0, (— /i ± \/h'^ — A2)^ — m^). These zeros fci,2,3,4 may 
not exist at the same time, but it can be shown that at least one of them must exist 
if there is a population asymmetry. As a consequence, a meaningful solution should 
have \h\ > A. It is also possible that the energy spectrum of fermions and that of 
antifermions may have different numbers of zeros. We also observe that the number 
of zeros in the energy spectrum of either fermions or antifermions may change. It is 
interesting that the topology of gapless regimes changes with the number of zeros in 
the spectra. However, this does not introduce any new phase or singularity. 

Number density of each species is determined by = \ Tr[7°G<7 {K)]. After 
summing fermion Matsubara frequencies and making use of Tr(7° A-t (k)7° ) = 2, one 
has 

= E + vi-^'fi-E-^) - {ui+^'mj + vi+^'n-Et,))] , 

k 

«i = E [<~^"m,) + vi-^'n-E-j - {ui^^'fiEi^) + vi-'^'fi-E^j)] , 

k 

(15) 



^ka 



A2, E+ = 



and for antifermions. 



and V 



(±)2 



-). Unlike 



from which on can determine contributions from fermions and antifermions. Here 
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f{x) = [exp(j:/T) + 1] is the Fermi distribution function. The polarization p is defined 

asp= (n| -nx)/(n| +ni). 

The F functions can be expressed as 



F^l{K,fi^,^fil) = TT^j 

,A-(k)^ 1 1 



75. (16) 



A similar expression can be found for Fii{K, From the definition of gap 

function, one can derive the gap equation. 

^ = -^ffETi-[75Fu(X,-/iT,Mi)] = -^.9ETi-[75nT(^,-A^i>m)]- (17) 



K K 



Explicitly, 



9 2Ei: 2E+ 



(18) 



We present another equivalent derivation of the gap function. There are two pos- 
sibilities for the asymmetric pair susceptibility in the presence of a number density 
asymmetry. They are defined as 

XTi(Q) = ^ETr[GT(^,A^T,-A^i)Goi(Q-if,m)] (19) 

K 

and xn (Q) with a similar expression. When Q = 0, it can be shown that 

XTi(O) = ^^Tr[»75Ftx(if,/iT,-/ii)], (20) 

K 

and Xn (0) has a similar expression. Furthermore, one already saw that these two 
expressions at Q = coincide. The gap equation follows Thouless criterion |[36|] which 
is given by 1 + gxix (0) = or 1 + gxn (0) = 0- They are equivalent to Eq. ( fTSl ). 

Since our simple model is not renormalizable, a proper regularization is required. 
Following Ref. LI 4.1 . the coupling constant is regularized by 



1ml X - / 1 



1 



g 47ra 2 V ej. — m + m 



(21) 



where a is the s-wave scattering length. We introduce a momentum cutoff A because 
there is an ultraviolet divergence. 

The nonrelativistic limit is defined as the limit when kpa- <C m, j/i^ — w| <C 
m, and A <C m. In this limit, all terms related to antifermion spectra or negative- 
energy solutions in the fermion propagator can be neglected. The energy spectrum of 
fermions becomes S^ka — ^ — {fJ-a — m) and the fermion propagator is reduced to the 
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usual non-relativistic results (for example, see OTIl ). Number densities also recover the 
nonrelativisitic results 13711 if contributions from antifermions are dropped. 
The ground state energy of the system is given by 

^ i^Tr{[7%A+(k)-7V.A_(k) + iET(^,-w)]GT(i^,m,-Mi) 

K 

+ [7°efcA+(k) -7"e,A_(k) + -^^t)] Mi, -Mt)} 

= E [^ka - + EUiEL) + EtpfiElp)] +^Jin + hSn - ^ . (22) 

k,A=± ^ 

The electrostatic energy is not included in our consideration. From the ground state 
energy one can derive the thermodynamical potential Qs = Es — TS — ii^n-\ — Mi'^i' 
which is 



f7c 
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E [ila-Et-TMl 



e T 



Tln(l + e-^)]. (23) 



k,A=± 

It can also be shown that 

n = — 



dfl J A,h ' 

give the same number equations and 



(^) 
\ dh J A, 



(^) =0 



(24) 



(25) 



gives the same gap equation. Our expression for is consistent with the result in 
Ref. ll28ll although in that work a variational ansatz with a coherent state is implemented 
to obtain fig. 



2.2. Bogoliubov Canonical Transformation 

The results obtained so far can also be derived from a canonical transformation. 
In previous subsections, our derivation follows conventional BCS ansatz, in which the 
ground state is assumed to be the BCS-Leggett ground state. Since the energy spectrum 
of quasi-fermions may become gapless, quasi-fermions survive in those regimes at zero 
temperature. In the deep BEC side when attraction is strong, the energy spectrum of 
quasi-antifermions can also be gapless. Thus in those regime quasi-antifermions can be 
found in the ground state. This is a new phenomenon which cannot be observed in non- 
relativistic systems. It is presumed that in BCS theory that fluctuations of the quantum 
operator gipJCi'^^tpi around its mean field value are negligible. As a consequence, the 
Hamiltonian derived from Eq. ^ can be linearized. It becomes 

H{x) ^ E • V + m - ^ct7°)?/;<t 

-A*(^f «75V'i) - A{^^zj,C^J) - — . (26) 



7 



The order parameter A can be taken as real with no loss of generality. Introducing 
annihilation and creation operators for fermions and antifermions, a„ (k), a\ (k), (k) 
and 6^(k), the fermion field operator can be expanded in the momentum space as 



1 



+ 



/2efe 



Ocr(k)e' 



-ik-x 



ik-x 



(27) 



where r}^ and r]^ correspond to two eigenvectors of a • k with different handedness 
and A;^ = (efe,k). 



cos^ 



sinacg^^k 



(28) 



Here and (}> are polar and azimuthal angels of the vector k. One can write the 
Hamiltonian in terms of a's and 6's. To make the result more compact, we redefine 
a<T(k) = e'~2~aCT(k) and 6cr(k) = e'~6o-(k), the Hamiltonian finally becomes 



^ k cr=T,i 

-A[ai(-k)aT(k) +6|(k)F4(-k) +4(k)fl|(k) + &i(-k)5T(k)]}. (29) 



Now we introduce the BCS-Leggett ground state for relativistic systems 

|BCS) = + v--a\{\,)-a\{-\,)][u+ + t;+6|(k)6| (-k)]|0) 



(30) 



and a BogoUubov canonical transformation 

/"(k)\ _ 



a(k) 
;3t(k) 



aT(k) 
at(-k) 

^T(k) 

t 

i 



(31) 



It can be shown that a(k)|BCS) = and /3(-k)|BCS) = 0. Therefore a(k) and 
/?(— k) are annihilation operators of a quasi-fermion and a quasi-antifermion, respec- 
tively. The Hamiltonian can be diagonalized following this transformation. It becomes 



H 



^ k 

+£;,-^/3t(k)/3(k) + £+ at(k)^(k) + i;+ ^t(k)^(k)) . (32) 
The order parameter in terms of this new set of operators has the following form 

^ = -9j2(^d-k)a^{k) + 6|(k)6t(_k)). (33) 
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With number densities for quasi-fermions (Q!'l'(k)a(k)) = /(E'^,^) and (/3'l'(k)/3(k)) = 
fi-^kp) those for antiquasifermions, it can be shown that Eq. ( |33] ) is equivalent to 
Eq. ( fTSl l. Similarly, number densities = (V'J-V'cr) give the same number equations 
as Eq. ( fTSl l. From the diagonalized Hamiltonian, one can calculate the thermodynam- 
ical potential Vis — —(3"^ liiTr (e^''^) , which has exactly the same expression as 
Eq. ( l23T l. Therefore these two formalisms are equivalent. We would like to point out 
that in the presence of inhomogeneity such as an external potential, the Bogoliubov 
canonical transformation can be generalized via the Bogoliubov - de Gennes equation 
1 38] to study these inhomogeneous phases. This will be particularly useful in the study 
of modulating phases such as the LOFF state. 



3. Stability Condition and Phase Separation 



After solving the number and gap equations consistently, one obtains a homoge- 
neous phase. However, it may not be stable. To check its stability, we analyze the 
number susceptibility matrix dna/dfia', which is related to the compressibility. The 
number susceptibility matrix reveals responses of the system to small perturbations 
in densities. Only when all eigenvalues of the number susceptibility matrix are non- 
negative, the phase is stable against phase separation in real space | Sfl 3^- This is the 
mechanical stability criterion. For a charged Fermi gas, there is another stability cri- 
terion which requires that the sqaured Meissner mass of the gauge bosons be positive. 
This additional stability criterion is discussed in BBIl for gapless color superconducting 
states. We only present results from the mechanical stability criterion. 

The gap function A is an implicit function of ii^r' (or ji and h) if the gap equation is 
satisfied. For simplicity in notation, we use the ordinary partial derivative ^ when A, 
/i and h are all treated as independent variables, and use -pj = ^ + ^ {x — /i, /i) 
as the total derivative with respect to /i and h when A can be inferred from the gap 
equation. The number susceptibility matrix can be written as 



Djii Dn 



/ Dn Dn \ 



-A 



D5n ^ 



Da 

Ddn 



1 



-AM A, 



(34) 



where A = A-^ = ^ 



Therefore, eigenvalues of the matrix are given 



by A± = Tr(Af) ± y/TiiM)^ - 4det(Af) /2, where Tr(Af ) and det(Af ) denote the 
trace and determinant of M, respectively. Matrix elements of M are list in AppendixiBl 
If a homogeneous solution is not stable, it is possible that phase separation in real 
space or a LOFF state may take place. It is known that the LOFF state is not likely to 
exist in the strongly attractive regime for non-relativistic fermions |27l] . Here we only 
consider phase separation for neutral relativistic fermions. In a phase separated state, 
there is an interface separating a superfluid and a normal phase. In the following sec- 
tions we will show that homogeneous superfluid phase is unstable in two regimes, one 
in the weakly attractive regime and one deep in the RBEC regime. Chiral condensates 
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with 7^ which are not considered in our model may also occupy those unstable 
regimes. 

Following ll40i 31], we consider a phase separated state with volume fractions of 



normal and superfluid phases being x and (1 — x), respectively. To achieve a stable 
phase separation, temperature, chemical potentials of the two species, and pressure 
should be continuous across the interface to balance heat convection, diffusion, and 
mechanical movement ||31[1 . Explicitly, we require Q.m —^s and = /i^. Here the 
thermodynamical potential of a normal phase is given by 

^N = -T [ln(l + e-^)+ln(l + e-^)]. (35) 

k,A=± 

When electrostatic energy is taken into account, a phase separated state may have 
higher energy due to charge accumulation. This may disfavor the phase separated 
state and other inhomogeneous phases such as the LOFF phase. In that case, other 
homogeneous phases (such as chiral condensates) should be considered. 



4. Results and Discussions 



Following 1 14], a BCS-BEC-RBEC crossover can be characterized by three dimen- 
sionless constants, the dimensionless coupling constant 77 — l/{kpa), ( = kp/m, and 
the cutoff z = A/{(m). The quantity C, describes how relativistic the system is. If 
C ^ 1, one should recover nonrelativistic results when the coupling constant is small. 
In this paper, we set A/m = 10 and focus on two cases: C = 0.1 (weakly relativistic) 
and C = 0.6 (relativistic). When 77 is extremely large, it is possible that p > 1 114 ill as 
one can see from Eq. (flSl l. Throughout this paper we consider p < 1. 

As mentioned in Sec. [T] relativistic effects become significant when rjc ^ V^C^^ 
for a single-species system. Close to this critical value of coupling constant, rjc, a 
RBEC-NBEC crossover is expected. We observe that r/c depends on p in the pres- 
ence of population imbalance. At zero temperature, we derive an analytical expression 
for the critical coupling constant. At rjc , zeros of energy spectra of fermions and an- 
tifermions can be found at /ci = = Q and k2 — k^ ~ \/h^ — vrp- by using the 
approximations /i ~ and A ^ m. It can also be shown that p ~ (kijkpY . With 
these expressions, the critical coupling constant is found to be 



2 r 



A 

771 



._A.2 



In I 



pa 



Cpi + 1 



(36) 



From this expression, one sees that when C ^ 1 in the nonrelativistic regime the first 
term dominates and therefore r\c is insensitive to p. For larger C, the second term is no 
longer negligible and r\c will have observable dependence on p. 

We demonstrate typical behavior of BCS-BEC-RBEC crossover in Fig.[T](for = 
0.1) and Fig. |2] (for Q — 0.6) by showing A, (i, and /i as a function of rj for selected 



10 










'(b) ■ 






























^ / 










0.4=1- / _ 





10 


20 / 


















0.95 








0.7 








— 0.05 



10 



15 



20 



Figure 1: (Color online) Signatures of a NBEC-RBEC crossover for f = 0.1. (a) A as a function of r] for 
selected values of polarization p (labeled next to each curve). Inset: A in the weakly attractive regime, (b) h 
and ^ (inset) as a function of r) for selected values of polarization p (labeled next to each curve). Solid dots 
indicate positions where relativistic effects become significant and crossovers take place. 




Figure 2: (Color online) Signatures of a NBEC-RBEC crossover for C, = 0.6. (a) A as a function of 7] 
for selected values of polarization p (labeled next to each curve), (b) h and fi (inset) as a function of r] 
for selected values of p (labeled next to each curve). Solid dots indicate positions where relativistic effects 
become significant and crossovers take place. 
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Figure 3: (Color online) Phase diagrams at zero temperature for (a) ^ = 0.1 and (b) f = 0.6. Phase transi- 
tions occur only when one crosses solid lines. NSF and RSF denote nonrelativistic and relativistic superfluid 
phases, respectively. Shaded regimes correspond to unstable superfluids, which are nom'elativistic in regime 
I and are relativistic in regime II. Black sold lines indicate where A = 0, which are boundaries between 
normal phases (N) and USFs. Blue dashed lines indicate rjc, which is given by Eq. ^36) and corresponds to 
the onset of RSF. Red dashed lines indicate boundaries between stable superfluids and unstable superfluids. 
Green dot-dash lines and maroon dot-dash lines indicate where the topology of gapless regimes in energy 
spectra of fermions and antifermions changes, respectively. To the right of the maroon dot-dash lines, gapless 
excitations of antifermions can be observed. 

values of p. Solid dots on those curves indicate where relativistic effects become sig- 
nificant. Their locations are obtained from Eq. (|36] l. In both cases, A, and h change 
dramatically near rjc- A and h become of order of Ep while /i approaches zero (and 
fermions and antifermions become nearly degenerate). One can see that locations of 
those solid dots shift to smaller coupling constants when the system is more relativis- 
tic. In both figures A decreases as p increases, this is similar to the nonrelativistic 
case since the presence of number density asymmetry frustrate pairing. We also found 
that /i has weak dependence on p while h has more observable dependence on p when 
77 > rjc- The reason for the former is that the system has a large two-body binding 
energy when the coupling constant is large. Since the chemical potential is dominated 
by this binding energy which is independent of p, the energy needed to break a bound 
pair (~ 2|/j, — m|) has a very weak dependence on p. The reason for the latter is that 
when the coupling constant is large, excess quasi-fermions form a sphere in momen- 
tum space similar to a Fermi sphere of non-interacting fermions. h is related to the 
radius of this sphere and therefore it is sensitive to Sn, which is controlled by p. These 
features are more prominent in Fig.[T] 

Fig. [3] shows a predictive phase diagrams in p-rj plane for (a) ^ = 0.1 and (b) 
C, = 0.6. When ( is small (weakly relativistic cases such as (a)), the phase diagram 
is similar to that of nonrelativistic cases, especially when 77 is small. We show loca- 
tions of on our phase diagrams. We found that crossing rjc does not introduce any 
phase transition. Therefore a nonrelativistic superfluid (NSF) and a relativistic super- 
fluid (RSF) are connected smoothly without singularities and we use dashed lines to 
indicate that only smooth crossovers occur. The dot-dash lines show where the topol- 
ogy of regimes of gapless quasi-fermions (green line) and quasi-antifermions (maroon 
line) changes. However, no new phase or singularities occur as the topology changes. 
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Figure 4: (Color online) Number density distributions (in unit of fc|,) for fermions as a function of momen- 
tum for p = 0.2 and C, = 0.1 with (a) rj = 19.8 (in the RBEC regime) and (b) t) = 10. Insets show 
corresponding number density distributions for antifermions. We use black and red to represent the two 
species. 




Figure 5: (Color online) Number density distributions (in unit of fcj.) for fermions as a function of mo- 
mentum for p = 0.2 and C, = 0.6 with (a) = 4 (in the RBEC regime) and (b) r} = —I. Insets show 
con'esponding number density distributions for antifermions. We use black and red to represent the two 
species. 



13 




ri ri 

Figure 6: (Color online) Phase diagrams including phase sepai'ation for (a) = 0.1 and (b) f = 0.6. Shaded 
regimes labeled by PS correspond to phase separation. For phase separation covering the unstable regime I, 
black solid lines indicate boundaries between normal (N) phase and phase separation and red sold lines indi- 
cate boundaries between stable nonrelativistic superfluid (NSF) and phase separation. For phase separation 
covering unstable regime II, violet solid lines indicate boundaries between stable relativistic superfluid (RSF) 
and phase separation. Brown dot-dash lines indicate where A = in unstable nonrelativistic superfluids. 
Blue and green dot-dash lines indicate where homogeneous superfluids become unstable. Orange dashed 
lines indicate rjc, which conespond to the onset of RSF. 



We also show unstable regimes of homogeneous superfluids and our results agree with 
those presented in Ref. |jl60. We also test the stability of homogeneous superfluids. 
As shown in Fig. Ob), there are two separated unstable regimes where number sus- 
ceptibility matrix is not positive-semi-definite. In addition to being unstable when the 
coupling constant is small, a relativistic superfluid is also unstable when the coupling 
constant is very large. 

In Fig.[3ja), the other unstable region of RSF appears when 77 > 59.7 and is not 
shown. We name the unstable regime with small values of r/ as regime I and that 
with larger values of ry as regime II. Inside the regime I, there are gapless superfluid 
states with two zeros in the energy spectrum of quasifermions, which are similar to 
the gapless CFL or 2SC states [42, 33]. The unstable regime II shifts deep into the 
RBEC side (with a larger coupling constant) as ( decreases (or the system becomes 
more nonrelativistic). Therefore observations of this unstable regime II is unlikely 
in nonrelativistic systems. As ( increases, the regime occupied by a nonrelativistic 
superfluid phase (NSF) shrinks. We found that it will eventually disappear completely 
at large enough ( and only RSF phases survive. In charged Fermi gases, it is likely 
that inside regime I, the Meissner mass of gauge bosons may become imaginary, which 
indicates a new unstable regime ||43ll . 

We present number density distributions as a function of momentum for fermions 
and antifermions with p = 0.2 in Fig. |4]for C — 0.1 and in Fig. |5]for ( = 0.6. 
When f] > rjc, antifermions develop a gapless regime in number density distributions. 
One can see that both fermions and antifermions are polarized. However, we observe 
an interesting feature: antifermions with spin-down has a larger population than an- 
tifermions with spin-up do, which is exactly opposite to the case of fermions. Due 
to this behavior of antifermions, the polarization of fermions needs to exceed the to- 
tal polarization of the system. This is another distinct feature of relativistic fermionic 
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systems. When the coupling constant is small, there is virtually no gapless excita- 
tion of antifermions. Comparisons between Fig.|4jb) and Fig.|5jb) show that the gap- 
less regimes of quasi-fermions have different topologies since the gapless regime in 
Fig Etb) has only one zero and the one in Fig.|5lb) has two zeros. 

In Fig.|6]we present predictive phase diagrams including phase separation (PS) for 
(a) C = 0.1 and (b) ( = 0.6. In previous sections, we found two regimes where homo- 
geneous superfluid phases become unstable. We call them region 1 with small values 
of rj and region 11 with large values of rj. In regime 1, LOFF states or phase separation 
may appear and compete with each other. In region 11, phases with chiral condensates 
(i^'ip) 7^ 0) or phase separation may appear In Fig. |6ta) only one of the two phase 
separation regimes is shown because the other one is too deep in the strongly attractive 
regime. For phase separation in regime 1, phase separation covers roughly the regime 
where homogeneous superfluids become unstable. In contrast, in regime 11 phase sep- 
aration covers a significantly larger regime than the regime where homogeneous su- 
perfluids become unstable. We also found that in regime 11, the volume fraction of the 
normal phase of a phase separation only increases slowly with the couping constant 
rj. For more complicated models with electrostatic energy, the boundaries of regimes 
where phase separation exists may shift. 

5. Conclusions 

In this paper we address relativistic BCS-BEC-RBEC crossover of fermions with 
tunable attractive interactions at zero temperature. We implement a BCS-Leggett ansatz 
and use two equivalent formalisms to investigate possible phases. We show that gapless 
excitations of quasi-fermions and quasi-antifermions are both possible. The stability 
of possible homogeneous polarized superfluid phases is analyzed. When homogeneous 
superfluid phases become unstable, we consider phase separation in real space and 
leave more exotic phases for future study. Number density distributions and predictive 
phase diagrams are presented to help understand effects of number density asymme- 
try on different phases in a relativistic BCS-BEC-RBEC crossover. Our results are 
consistent with the well-known results in the nonrelativistic Umit. 

It is of interest to extend our study to finite temperatures. It is conjectured that color 
superconductivity vanishes at high temperatures and the system becomes a quark gluon 
plasma. From study of BCS-BEC crossover, it is known that condensed Cooper pairs 
can either be broken into quasi-particles or be excited and become non-condensed pairs 
[44]. To include both possibilities, one needs to consider pairing fluctuations and this 
is beyond mean field theories. In Ref. ifisll finite-temperature effects are considered 
for a one-species Fermi gas. It is possible to extend this approach to include finite- 
temperature effects for a two-species relativistic Fermi gas in future work. 

This work is supported by Grant NSF PHY-0555325. We thank K. Levin for helpful 
discussions and P. Zhuang and L. He for useful comments. 
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A. Differential operator matrices 

The differential operator matrix is 



D{x) = 



(37) 



iA75 

iA75 Doi{x,-^i{) 

Here Doa{x, i^a ) = — 7°i9t + ^7 • V — (to — /Xcr7'') is the noninteracting differential 
operator matrix. In momentum space, 

G-/(is:,MT) iA75 



D{K) 



iA75 G^l{K,-ni) 



(38) 



where Go^^-ftT, /x^) 
propagator. 



{iuJn + Mo-)7 — 7 • k — m is the inverse of a non-interacting 



B. Number susceptibility matrix 

The elements of the transformed number susceptibility matrix are given by the 
second-order derivatives of the thermodynamic potential 

Dn__dHls /a^y 
Dn _ d'^Qs d'^ils d'^^s 5^0.<? DSn 



DSn 



d^idh 



s 



dudAdAdh' aA2 Dji 



(39a) 
(39b) 



Dh dh^ \dAdhJ ' dA^ ' 
The expressions of the related second-order derivatives are 



(39c) 



dA^ 

dHls_ 

dAdu 
dAdh 
dfidh 

F+ = 



EZ 



A2 A-p 



k 



/- 



E. 



Vet^ ^ e+^ 



+ /'-] + 



Ef 



Et 



Ef 



Ucrcf+^fiE^J+fiE^^), r 
f-^fiE^J+fiE^0),f'+-. 
and f'{x) = df(x)/dx. 



^f{E^o.)+f{Eu0ir 

f'iEtj-nE+f,),!' 



f'iE^J 



(40) 

fiEt^), 
fiE^f,), 
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